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We study the influence of a strong imaginary vector potential on the quantum mechanics of 
particles confined to a two-dimensional plane and propagating in a random impurity potential. We 
show that the wavefunctions of the non-Hermitian operator can be obtained as the solution to a 
two-dimensional Dirac equation in the presence of a random gauge field. Consequences for the 
localization properties and the critical nature of the states are discussed. 
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The field of non-Hermitian quantum mechanics has at- 
tracted great interest recently both in its connection to 
anomalous diffusion in random media |^ , as well as to the 
statistical mechanics of flux lines in superconductors . 
At the same time, motivated in part by the connec- 
tion to properties of the integer quantum Hall transition 
and gapless superconductors detailed investigations have 
been made into the critical properties of Dirac Fermions 
coupled to random gauge fields [p[-pl|. The aim of this 
letter is to identify a connection between these prob- 
lems which explains some of the unusual phenomena re- 
cently reported in the behavior of two-dimensional non- 
Hermitian random operators This correspondence 
is related to a chiral symmetry of effective Hamiltoni- 
ans commonly used in the analysis of problems in non- 
Hermitian quantum mechanics. 

The two-dimensional Hamiltonian we consider de- 
scribes a particle propagating in a random scalar im- 
purity potential, V and subject to a uniform imaginary 
vector potential, ih 
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The scalar potential, V is assumed to be real and drawn 
from some random distribution, P[V] which, for now, is 
left unspecified. Related problems have been recorded in 
a variety of physical situations ranging from the study of 
reaction diffusion phenomena in biological systems [p^ , 
to advective diffusion in random media |^Jl^,|lJ] , and the 
study of fluctuating vortex lines in superconductors with 
columnar defects 

Remarkably, in contrast to properties of the Hermitian 
operator (i.e. one in which the vector potential is real), 
numerical studies |^ suggest localization properties of H 
depend sensitively on the relative strength of h. While, 
in dimensions d < 2, all wavefunctions of the Hermitian 
operator are believed to be localized [ p^ , numerical ev- 
idence suggests that application of a sufficiently strong 
vector potential, h induces a delocalization transition of 
states of the non-Hermitian operator. In contrast to the 



situation in Id PlJlq] , the mechanism and stability of the 
delocalization transition is not yet well understood. 

Recent analytical studies have focussed on the spec- 
tral properties of the non-Hermitian Hamiltonian [ pT| . 
Treating the vector potential as a weak perturbation of 
the random Hamiltonian, impurity averaged properties 
of the Green function have been cast in the form of a 
functional field integral involving a supcrsymmctric non- 
linear CT-model. While capturing universal features of 
the complex spectrum of H, evidence for delocalization 
of states was not sought. By contrast, in the present 
approach, we will impose a strong imaginary vector po- 
tential, and treat the random potential as a perturbation 
(a regime explored in Id by Feinberg and Zee |ll6|] ). By 
doing so we will reveal an explicit connection between 
the Hamiltonian in Eq. (^ and the problem of Dirac 
fermions propagating in a random gauge field. In par- 
ticular, we will find a regime of weak disorder in which 
it is possible to construct explicitly eigenfunctions of the 
non-Hermitian Hamiltonian H for individual realizations 
of the disorder. 

In the absence of the impurity potential, the eigen- 
functions of H are plane waves with complex eigenval- 
ues, zo{p) — (p^ — h^)/2TO + ip ■ h/m. Since the Hamil- 
tonian is real {H = H*), eigenvalues occur in complex 
conjugate pairs, a property maintained in the presence 
of the random potential. In the infinite system, the spec- 
trum forms a dense support occupying the region of the 
complex plane |Im z\ < y/2mRe z + |hp |h|/m. The 
two-dimensional density of states (DoS) takes the form 
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[(2mRe z h2)h2 - TO2(Im z)^]^/^ ' 
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where denotes the constant DoS of the Hermitian 
Hamiltonian. 

Concerned with impurity averaged spectral properties 
of the Hamiltonian, we begin by defining the Green func- 
tion G{z) = {z — H)~^ where z denotes the complex 
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argument. Using = (l/7rri)tr dz'G{z), where f2 rep- 
resents the volume, G{z) is shown to be non-analytic ev- 
erywhere the DoS z^(z) is non- vanishing Jl^ , p^ . 

To properly account for non-analytic properties of 
the im purity averaged Green function, previous stud- 
ies have emphasized the need to express the 
Green function through a Hamiltonian which is explic- 
itly Hermitian. This is achieved by constructing a matrix 
Hamiltonian with the 2x2 block structure 



n = 





iJt - z* 



H - z 
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In this representation, the Green function of the non- 
Hermitian operator is expressed as the off-diagonal 
element of the matrix Green function, G{z) = 
lim^^o G2i{z), where, defining ij ^ 0+ ^ Q — {irj — 'H)~^ . 
Zero energy eigenstates of the matrix Hamiltonian Ti. 
yield eigenstates of the non-Hermitian Hamiltonian H. 

Defining Pauli matrices a ((Tq = 1) which operate in 
the 2x2 space, the chiral symmetry of the matrix Hamil- 
tonian, Id — —a^Tia^ implies that eigenvalues, of Ti. 
appear in pairs of opposite sign. Moreover, any such pair 
of eigenstates obeys \ip+ei) — o'slV'-ci)- 

The spectrum of Ti. depends sensitively on the strength 
of the vector potential, h. In the absence of the random 
potential, the dispersion relation takes the form. 



Eoip) = ±\zo{p) 
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invariant under reflection about the axis parallel to h. 
Thus, in contrast to a Hamiltonian involving a real vec- 
tor potential, where the spectrum is described by a simple 
shift of the Fermi sphere, the continuous degeneracy of 
the zero eigenvalues (the poles of the Green function) is 
lifted. Instead, setting h = he2, zero energy states exist 
at only two discrete points Pq''^ a — 1,2 (see Fig. |^). 

If the impurity potential is strong (i.e. i <C h/h, where 
£ denotes the transport mean free path associated with 
the random impurity potential, V), unperturbed states of 
the clean system are strongly mixed by the disorder (see 
Fig. In this limit, the pole structure of the impurity 
averaged Green function, Q is smeared out. Correspond- 
ingly, statistical properties of Ti. are largely insensitive 
to the zero eigenvalues of i?o(p). In this limit, one can 
expect the transport properties of the non-Hermitian op- 
erator to reflect those of the Hermitian counterpart. Con- 
versely, if the impurity potential is weak {i 3> h/h) the 
pole structure of the average Green function is dominated 
by the nature of the spectrum in the vicinity of the zeros, 

(a) 

Po ■ 




FIG. 1. Dispersion relation, -Eo(p) (normalized by vih) in 
the vicinity of zero energy shown as a function of p with 
h = he2, Im 2 = and mvi/2h = 10. Note that the mini- 
mum scattering amplitude required to smear out the points 
of degeneracy is given by 7i/r ~ hvi. 



Focusing on the limit of weak disorder 

Re 2 > - — > - , Im z, 
Zm T 
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where r represents the corresponding mean free scatter- 
ing time of the random potential, a linearisation of the 
spectrum in the vicinity of the two zero eigenvalues sepa- 
rates the spectrum into two branches. Treating the ran- 
dom potential and Im z as weak perturbations, and per- 
forming a gradient expansion in 
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where Pq"'' = {—l)'^^^mvie^, mvi = (2TORe z + ft,^)^/^, 
and mv2 = h, the low energy, long-wavelength expan- 
sion of the (unperturbed) Hamiltonian around the Fermi 
points Pg°'* generates the anisotropic Dirac operator 
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-To (Xl [criWi(p • ei) + (J2V2{p ■ 62)] 
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Here we have introduced an additional set of Pauli matri- 
ces, r (tq = 1), that index the block structure associated 
with the reflection symmetry. The existence of two de- 
generate zero energy eigenstates of 7i^^ conspires with 
the exact chiral symmetry to yield an anisotropic Lorentz 
symmetry. 

Being generally non-symmetric under reflection, ma- 
trix elements of the random impurity potential violate 
the reflection symmetry. Accounting for matrix elements 
which scatter across the Fermi surface (i.e. between the 
Dirac points) as well as within each subspace, the general 
Hamiltonian takes the form 
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Ud^T-(-d^ +^'ry®<7iyiy+To(^a2lTa z, (8) 

where the impurity potentials, Vi/(r) are real random 
functions with Fourier components 



V{-qi,q2) V{qi~2mvi,q2) 
V*{qi - 2mvi,q2) V{+qi,q2) 
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Equations represent an important intermediate re- 
sult of this Letter: Firstly, the low energy sector of the 
original Hermitian Hamiltonian has been described 
in terms of the stochastic Dirac Hamiltonian, TLd- Sec- 
ondly, ItLd possesses the chiral symmetry cr^HDCa = 
—Til), a direct consequence of the chirality of the auxil- 
iary operator (||). The significance of the second obser- 
vation lies in the fact that the behavior of chiral Dirac 
Hamiltonians can be analyzed for individual realizations 
of the disorder. As we shall see, it is possible to construct 
explicit solutions for the zero energy eigenfunction I'^'d) 
of Hd- From these the eigenfunctions jvPz ) (j^'x.)) of 
the original non-Hermitian operator H (i?^) can then be 
obtained as 







ffl'l*, 







D 
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for an arbitrary eigenvalue z provided Eq. (^) holds. Here 
the matrix structure refers to the a-space and the opera- 
tor H accounts for the fact that in order to obtain eigen- 
functions of the full problem, the eigenstates of the low 
energy expansion (^) have to be 'boosted' to the Fermi 

points Po"\ a — 1, 2, respectively 

We next turn to the explicit construction of the 
wavefunction ^!]j{r). To this end we first remove the 
anisotropy of the Hamiltonian by rescaling the coor- 
dinates according to = {v^/v)x^, v —^v\V2- As a 
result Tiu takes the canonical form 



'Hd 
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where h = v = 1 and the two components (/i = 1,2), 
Afj. ^ + Bf, + are given by 
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Af_, = df_,iToVa, B^ = 5^1 ^ r^^K, C"^ = '5;i2ToIm z. 



The disorder in Ti/j appears in the form of a minimal 
coupling to a generally non-abelian vector potential (i.e. 
B^ 7^ 0). It is thus natural to seek a gauge transfor- 
mation that removes the stochastic components of the 
Hamiltonian. Indeed, although the potentials, A^j_ and 
B^ are not in general of pure gauge type, the non-gauge 
components can be accounted for by extending the con- 
cept of gauge transformations so as to include 'axial' 
transformations. 

Focusing on the abelian sector first, we decompose A^ 
into a transverse (axial gauge) and a longitudinal (pure 
gauge) component: A^ = e^^di,x± + d^x\\i respectively 
Ipif . It is then straightforward to verify that 



*d(x) 



_ iCToa:2lmx iCToXll(x) (T3Xi(x) 



e+ 
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represents a solution oilhio'^D = for = (22|. Here, 
Q± € C and the gauge transformation exp(i(To2;2lmz) 
has been used to dispose of the small imaginary compo- 
nent of the eigenvalue z in the non-Hermitian problem. 
This contribution can be absorbed into the 'boost' com- 
ponent of H (|) m. 

The treatment of the non-abelian components B^ 
is conceptually similar but - due to their non- 
commutativity in r-space - technically more involved. 
Referring to Ref. for details, we merely state that in 
the presence of finite B^, Eq. ( p^ ) generalizes to 



*d(x) = 



f{xi + iX2) 6 + 

g{xi - ix2) 6- 



C/(x) e 



c'3X-L(x)„o-3C(x)-r 
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Here, ?7(x) g SU{2) plays the role of the abelian gauge 
transformation exp[icroX||(x)], whereas the vector field ^ 
plays the role of the abelian axial component x± 
obeys an analogous but more complicated equation [ p4| . 
Finally, / and g represent analytic functions which are 
fixed by the boundary conditions (see below). 

Equation ( pT| ) represents a non-perturbative solution 
to the zero energy Dirac equation for any given realiza- 
tion of the disorder. It, therefore, allows the construc- 
tion of the eigenfunctions of the non-Hermitian Hamil- 
tonian according to Eq. (^. The most important prop- 
erties of are a) that the generalized gauge factors 
depend in a non-local way on the spatial distribution of 
the disorder potential and b) that the axial gauge factors 
lead to an exponential disorder dependent amplification 
of the modulus of the wavefunction, ^E'l^'J'jj = I'I'jjP = 

9tg2CT3 [x±+i-r]Q^ 

What can be said about the asymptotic behavior of 
the wavefunction and, in particular, about its lo- 
calization properties? To address this question we con- 
sider the impurity averaged two-point correlation func- 
tion, C(ri - ra) = |*D(ri)p9i |*£,(r2)|2'?2. From the 
structure of the solution we infer that C — Fa x Fb 
factorizes into an abelian and a non-abelian component 
both of which can be straightforwardly extracted from 
Eq. (pT|). To say more, it is necessary to specify both the 
form of the random potential distribution function, and 
the topology of the system. Here we focus on the ther- 
modynamic limit with Gaussian white-noise distributed 
disorder, P[V] of uniform variance. Under these circum- 
stances, taking / = 5 = 1 in Eq. ( pT| ) represents the only 
admissible choice [ p5[ . 

Defining gA ~ (2 x 2ttvt)^^ as the variance of 
the coarse-grained Gaussian distribution for (i.e. 
Xj^(x)x_L(y) oc — (7Aln|x — y|), and focusing, for sim- 
plicity, on the case qi =92 = 1, the impurity average 
yields a strongly anisotropic (wi 3> V2) algebraic decay, 



Fa{v) o,\{virif + {V2r2f\- 



(12) 
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implying that eigenstates of the non-Hermitean operator 
H are critical Moreover, since the algebraic decay 

of Fa is not influenced by Im we infer that the critical 
nature of the wavefunctions is also insensitive to Im z. 

As for the non-abelian sector, the non-trivial relation- 
ship between the fields and the effective 'gauge' field, 
^ makes the calculation of the correlation function, Fb 
more involved (c.f. Refs. However, although 

at present no rigorous statements on the long distance 
behavior of the correlation function Fb can be made in 
general, insight can be drawn from the following facts: 
a) The theory possesses a strong coupling fixed point 
described by the Wess-Zumino-Novikov-Witten theory 
SU -4^{2) (corresponding to an infinite strength of 

the non-abelian components of the disorder potential), 
and b) in this limit [0,0 , the correlation function, 

FB{v)o,\{v^r,f + {v2r2f\-''^ (13) 

is again algebraic. Note that the scaling exponent in Eq. 
( |l3| ) is fixed solely by the number of nodes (two) in Eq. 
(4). Therefore, given that the wavefunctions are critical 
in the strongly disordered limit, it seems highly plausible 
that they remain critical in general. 

In this Letter, we have studied the spectral properties 
of the two-dimensional random Schrodinger operator in 
the presence of a uniform imaginary vector potential, h. 
Mapping it to a Hermitian Hamiltonian with chiral sym- 
metry, a gradient expansion identifies properties of the 
non-Hermitian operator with those of a stochastic Dirac 
Hamiltonian Ti d ■ This correspondence allows for the ex- 
plicit construction of eigenfunctions of the non-Hermitian 
Hamiltonian for individual realizations of the disorder. In 
the thermodynamic limit, the wavefunctions were shown 
to be delocalized along both the directions parallel and 
perpendicular to h. 

Finally we note that the main characteristics of the 
wavefunction (|l^), long-ranged disorder dependence en- 
coded in the 'gauge fields', and exponential amplification 
of the wavefunction modulus, can manifest themselves 
in more complex phenomena than that discussed in this 
Letter: Firstly, one can envisage systems with non-trivial 
topology and/or stochastic potentials with superimposed 
regular structures (e.g. spatially non-uniform distribu- 
tion functions). In such cases the behavior of the wave- 
functions may change qualitatively (e.g. to localization). 
Secondly, we note that the sensitivity of the modulus to 
disorder results in strong statistical fluctuations of the 
wavefunctions which are large in comparison to the aver- 
age. In particular, correlation functions of the moments, 
C acquire scaling exponents with non-linear dependence 
in qi, a characteristic related to multi-fractality. 

We are indebted to A. Andreev, P.W. Brouwer, J.T. 
Chalker, D.R. Nelson, D. Taras-Semchuk and X.-G. Wen 
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from the Swiss Nationalfonds. 
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